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Abstract 

TBA integral equations are proposed for 1-particle states in the sausage- and SS-models 
and their a-model hmits. Combined with the ground state TBA equations the exact mass 
gap is computed in the 0(3) and 0(4) nonhnear cr-model and the resuhs are compared 
to 3-loop perturbation theory and Monte Carlo data. 



1 Introduction 



A better theoretical understanding of finite size (FS) effects is one of the most important prob- 
lems in Quantum Field Theory (QFT). The study of FS effects is a useful method of analysing 
the structure of QFT models and it is an indispensable tool in the numerical simulation of 
lattice field theories. 

Liischer derived a general formula for the FS corrections to particle masses in the large 
volume limit. This formula, which is generally applicable for any QFT model in any dimension, 
expresses the FS mass corrections in terms of an integral containing the forward scattering 
amplitude analytically continued to unphysical (complex) energy. It is most useful in 1 + 1 
dimensional integrable models [21, where the scattering data are available explicitly. 

The usefulness of the study of the mass gap in finite volume is demonstrated jHj by the 
introduction of the Liischer- Weisz- Wolff running coupling that enables the interpolation be- 
tween the large volume (non-perturbative) and the small volume (perturbative) regions in both 
two-dimensional sigma models and QCD. 

An important tool in the study of two-dimensional integrable field theories is the Thermo- 
dynamic Bethe Ansatz (TBA). This thermodynamical method was initiated by Yang and Yang 
jl] and allows the calculation of the free energy of the particle system. The calculation was 
applied to the XXZ model by Takahashi and Suzuki j3j who derived the TBA integral equations 
for the free energy starting from the Bethe Ansatz solution of the system and using the "string 
hypothesis" describing the distribution of Bethe roots. 

The TBA equations also determine FS effects in relativistic (Euclidean) invariant two- 
dimensional field theory models where the free energy is related to the ground state energy 
in finite volume by a modular transformation interchanging spatial extension and (inverse) 
temperature. Zamolodchikov [Hj initiated the study of TBA equations for two-dimensional 
integrable models by pointing out that TBA equations can also be derived starting from the 
(dressed) Bethe Ansatz equations formulated directly in terms of the (infinite volume) scattering 
phase shifts of the particles. In this approach the FS dependence of the ground state energy 
has been studied fl^. in many integrable models, mainly those formulated as perturbations of 
minimal conformal models. 

The TBA description of excited states is less complete. The excited state TBA systems first 
studied jHl IH] are not describing particle states, they correspond to ground states in charged 
sectors of the model. An interesting suggestion is to obtain excited state TBA systems by 
analytically continuing [Hj those corresponding to the ground state energy. TBA equations for 
scattering states were suggested for perturbed field theory models by the analytic continuation 
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method jTU]- Excited state TBA equations were also suggested for scattering multi-particle 
states for the Sine-Gordon model at its = 2 supersymmetric point [TTj . 

In [12] we proposed TBA integral equations for the excited states in the Sine-Gordon (SG) 
model (and massive Thirring (MT) model). Although in the SG/MT case the excited state 
TBA description is "superfluous" since based on the Bethe Ansatz solution of the model we 
already have the Destri-deVega (DdV) nonlinear integral equations [13 CH CSj to study FS 
physics, the simple pattern of the excited state TBA systems we found there seems to suggest 
that similar systems can be found also for other models, where no DdV type alternative is 
available. 

Our aim in this paper is to calculate the finite volume mass gap in the 0(3) and 0(4) 
nonlinear a-models. These can be represented as some special limits of well-known integrable 
models: the sausage-model [TE] and the SS-model [TTj respectively. The ground state TBA 
equations are known for both models and although no Bethe Ansatz solution is available, 
based on our SG experience we make the following assumption: excited state TBA equations 
and Y-systems exist in these models and they are (almost) of the same form as for the ground 
state problem. 

To transform the Y-system equations into TBA integral equations we need to know the 
analytic properties of the Y-system functions, in particular the distribution of their zeroes. 
Using Liischer's asymptotic formula we can calculate this distribution in the infinite volume 
limit. Our second assumption in this paper is that the qualitative properties of this distribution 
remain the same for finite volume. Using this conjecture we can write down the complete set 
of TBA equations that are sufficient to calculate the finite volume mass gap for the sausage- 
and SS-models and their cr-model limits. 

If both assumptions are true, the solution of the TBA problem provides the exact value of 
the mass gap. We have numerically computed the mass gap for both the 0(3) and the 0(4) a- 
models and compared them to Monte Carlo (MC) results and 3-loop perturbation theory. The 
agreement with asymptotically free perturbation theory (PT) (for small volume) is especially 
important since our starting pont was Liischer's formula (valid for large volume). 

The paper is organized as follows. In section 2 we summarize the S-matrix data for the 
integrable models that are considered in this paper. In section 3 we recall the TBA integral 
equations and Y-systems corresponding to the ground state problem. In section 4 we briefiy 
summarize the results of [12] for the TBA description of excited states in the SG(MT) model. 
In section 5 we apply Liischer's asymptotic formula to the integrable models of this paper 
and this is used in section 6 to write down the full infinite-volume solution of the Y-systems 
of the sausage- and SS-models and their a-model limits. The complete TBA description of 
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the 1-particle states of the SS- and sausage-models are given in sections 7 and 8 respectively. 
Numerical solution of the TBA integral equations is discussed in section 9 and the results are 
compared to available MC and PT data in section 10. Finally our conslusions are summarized 
in section 11. 



2 S-matrix data 

In this section we briefly summarize the S-matrix data and some other properties of the models 
which will be considered in this paper. 

The SG model 

We first consider the Sine-Gordon model and parametrize the SG coupling as 

= (1) 

For p > 1 we are in the repulsive regime and a soliton |+, ^) and an antisoliton| — , 9) of mass M 
form the spectrum of the model. The bootstrap S-matrix of the model is as follows [T8] : 



sXX{e) = szzie) = Aie) = -explzxi-] \ , (2) 



TT 



00 



I dk . ,^ ^. sinhfw — 1)A; 

x(0 = 2 / -r MkO 9ik), g{k) = - — , .J (3) 

k 2 cosh(/ej ■ smh^p/c) 







StZ{9) = Szm = K B{e), (4) 



sinh - 

Bi9) = Ai^n-9) = -—^Aie), (5) 

p 

sinh 



27r 
P 



Sl-ie) = 3-1(9) = 0(6) = — Ai9), Citn -6) = 0(6). (6) 

p 

This S-matrix is a solution of the Yang-Baxter equation, has a U{1) symmetry and satisfies 
the usual requirements of C-, CPT- and Bose-symmetry, unitarity and real analiticity for k = 
±1. The sign difference between the two possible choices of k, becomes relevant in the crossing- 
relation: 

Sl'Mn-e) = Cp,Cs.SZi0), (7) 
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where Cap is real and symmetric for k = +1 and Cap is imaginary and anti-symmetric for 
K = —1. The physical SG model corresponds to the choice k = +1 but we note that the above 
S-matrix becomes S'f/(2)-symmetric in the limit p —>■ oo for the choice k = —1 . 



The SS-model 

The next model we consider is the SS-model ^j. There are four fundamental particles of mass 
M in its spectrum: \A,6) , where 

A={ai,a2) ai,a2 (8) 

The S-matrix of the fundamental particles J7] is: 

'§^fw = -^2:lw^2Sw, (9) 

where Slp{9) is the SG S-matrix with parameter p and S]^p{0) is the SG S-matrix with parameter 
p, with the same k parameter value. For p,p > 1 we are in the repulsive regime and there are 
no bound states of the fundamental particles in the model. This S-matrix satisfies the crossing 
relation 

S^Ei^^ -d) = CbmCdnS'1^\0), (10) 
where Cab is a real symmetric matrix with the following non zero matrix elements: 

= = c'(+,_)(_,+) = = 1. (11) 



The sausage-model 

There are 3 fundamental particles of mass M: \a, 9) , a E {+, — , 0} in the sausage-model jTHj. 
When the coupling constant < A < |, there are no bound states in the model. The S-matrix 
elements of the fundamental particles fH] are: 



S%iO)^St(^.-0)^-. J^^^^^^ .Stm, (13) 



q+-(0) = sin7rA-sin27rA 

+ ^ ^ sinhA(^-2z7r) ■sinhA(e + z7r)' ^ ' 
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(15) 



SM = Stm + S^-{9). (16) 
The 0{n) non-linear cr-model 

The 0{n) NLS model consists of n self-conjugate particles: \a,9) of mass M , a G {1,2, ...,n}. 
The S-matrix of the model is 



where 



SM = (T,{9)6ab6cd + Cr2{9)6acSM + Cr3{9)6adSbc, (17) 

2^9 s^^\9) 



a,{9) 



in -9 {n-2)9-2ni'' 



(19) 

aJ9) = -2m ■ (20) 

' {n- 2)9 -2m ^ ' 

and the 'isospin 2' phase shift s^^-* is given by 

s(2)(0) = -exp J 2« j ^sm{iu9)-kniuj) I (21) 







with 



-7TU1 I p-^TT^ 

^n(-) = ■ (22) 

The concrete values of the 'isospin 2' phase shift for some low values of n are: 

n = 2 s(2)(0) = -exp|^xoo(^)}, (23) 
n = 4 s(2)(0) = _exp|2^xoo(^^^|, (24) 

where Xoo(0 is the p ^ cx) limit of x(0- 
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These models can also be obtained from the models discussed previously by a limiting 
procedure. Concretely, the 0(2) model can be obtained from the SG model in the limit p —>■ oo 
|19| 1^ . the 0(3) model from the sausage-model in the limit A — ^ ^Hl and finally the 0(4) 
model from the SS-model (with k, = —1) in the limit p,p ^ oo In the rest of the paper 

only the k = —1 SS-model will be considered and called SS-model. 

3 Ground state TBA equations and Y-systems 




b 



Figure 1: Dynkin-diagrams associated with Dp+i and A2p_i type Y-systems. 



In this section we will give a short review of the ground state TBA equations and Y-systems for 
the models considered in the previous section. The TBA equations of these models (at some 
special values of the coupling constants) can be encoded in a 'Dynkin-diagram'. The unknown 
functions Ya{C,) are associated to nodes of the Dynkin-diagram and the TBA equations are of 
the form 

Ya{^) = e-''^"i^°'^'^f« e'^'"^^), / = ML, (26) 

where 

la{u) = J2lab\n[l + Y,{u)], (27) 



oo 



— oo 

M is the mass of the particles, L is the box size and lab is the incidence matrix^ of the Dynkin- 
diagram. 

The ground state energy can be calculated from the solutions of the TBA equations : 

oo 

^(o)__^ y rfu cosh In [1 + Fi (m)] . (29) 



^lab is zero if nodes a and b are not connected by links and it is unity if the nodes are connected by a single 
line. 
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p -1 p -2 2 1 2 p-2 



Figure 2: Dynkin-diagram associated with the SS-model Y-system (ground 
state) . 



The SG TBA equations for integer p >2 correspond to the -Dp+i type Dynkin-diagram shown 
in Figure la pi]. The TBA equations of the SS-model (for integer p,p>2) correspond to the 
Dynkin-diagram shown in Figure 2 ^^I, finally the TBA equations of the sausage-model (for 
A = with integer) are associated to the Dynkin-diagram shown in Figure 3. 
The solutions of these TBA equations are also solutions of the so called Y-systems 



(30) 



The standard way of solving the TBA equations is to iterate starting from the large / 
solution. The leading / oo solution of the TBA equations is easily obtained for the models 
discussed above. The solutions are listed below. 



SG model 



YAu) 



2g— ' cosh 



(31) 



Dp constant solution 



SS-model 



Yk{u) 
Yp{u) 



Y 



p+i 



u 



k = 2, 
p — 1 



,p-l 



YUu) ^4e 



'I cosh 



(32) 



(33) 



Sausage-model 



Yk{u) : Dp constant solution k = 2,...,p+l 



Yk{u) : Dp constant solution k = 2,...,p + l 



Yk(u) : Dn constant solution k = 2,...,N + l 



(34) 
(35) 
(36) 
(37) 
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Figure 3: Dynkin-diagram associated with the sausage-model Y-system. 



4 Excited states in the SG model 

In this section the excited states TBA equations ^2] of the SG model (and the closely related 
massive Thirring (MT) model) will be briefly summarized. The SG (MT) model can be regu- 
larized on a light-cone lattice in an integrable way The regularized lattice model can be 
solved by the Bethe Ansatz method. From the Bethe Ansatz solution of the model it follows 
that there exists a Y-system of the form of pO|l for all excited states of the model. For sim- 
plicity in the rest of the paper we consider only the p > 2 integer case with H soliton/fermion 
states without antiparticles. When H is even or both H and p are odd then the corresponding 
Dynkin-diagram is of -Dp+i type. In the latter case if p > if then the -Dp+i diagram is reduced 
(Yp = Yp^i = —1, Yp^i = 0) to the Ap^2 type diagram shown in Figure 4. When H is odd and 
p is even then lab is the incidence matrix of an v4|p_^ type diagram which is shown in Figure 16, 
where the oriented double line at the end of the diagram means 

Ip~ip = 1, Ipp-i = 2. (38) 

One can see that the same Y-system pup describes a large number of different excited 
states of the model. The difference between the various excited state solutions of the Y-system 
is in the analitical structure of the solutions. The Y-system functional relations (jHUj) can be 
translated into TBA integral equations in a standard way For this we have to know the 
positions of the zeroes and poles of Ya{C,) in the strip |lm,^| < 1 together with their asymptotic 
behaviour. We call this strip the main strip. From the Bethe Ansatz solution it follows that 
the Ya{C,)'s can have only zeroes in the main strip. 

The set of zeroes of Xi(0 (i^ the main strip) will be denoted by 

1a = H"^}ti- (39) 
These zeroes are related to the T-system zeroes 

-^ = W^^£- (40) 

as follows. 

gi = qa = ra-iUra+i a = 2, (41) 
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(a = 2, ... ,p — 3 for the Ap_2 case) and 

= 2-r,_i (42) 

qp^2 = rp^3 P > 5). 

With these definitions the TBA integral equations are of the form 

Qa 

a=l 



where cTq is the sign of Ya{oo) and r(0 = tanh (f ^). 



-o — o o — o 



1 2 3 p-3 p-2 

Figure 4: Dynkin-diagram associated with Ap^2 type Y-systems. 



These equations have to be supplemented by the "quantization conditions" 

1 + Ya (yi") ± z) = n = l,...,Ra a = l,...,p. (44) 

(a = 1, ... ,p — 2 for the Ap_2 system.) The exponential factor e"'™*^^!^ is present in the TBA 
equation for a = 1 only. This is indicated in the figures by colouring the corresponding nodes 
black. 

An important special case is when all zeroes are real. In this case the modulus of Ya {y^^ ± i 
is automatically equal to unity and ()44j) can be rewritten as 



(^)^" exp -I sinh - a, (yi")) + g 7 [v^a^ - 4"^) | = "^a, (45) 

for n = 1, . . . , i?a and a = 1, . . . ,p (a = 1, . . . ,p — 2 for Ap_2). In (pSj) the notation 

7(m) = 2 arctan (t(m)) , (46) 

00 

«,(«) = / rft; . , , (47) 

4 J smh ^(t> — 



-00 
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is used, where V indicates principal value integration. Note that |7('u)| < | for real u. The 
quantization condition for Yi{^) 







n 



Ri 



plays a special role. The set of solutions of (jlHjl is 



ri 



Ri 



n=l 



(4J 



(49) 



which includes the real zeroes {ha}a=i (Bethe Ansatz 'holes') and complex zeroes {^p}^^-^ 
(complex 'holes'). 

The energy and momentum of the model can be easily expressed by these zeroes and by 

H , C 



E 



M 



M 



cosh — ^ + cosh 

/3=1 



/3 



a=l 
■ H 



sinh 



2 

TChr 



a=l 



C 

+ sinh 

(3=1 



2 



TCU 



du cosh — ln[l + l^i('u)] 



1 /"°° TIU 

- I (iwsinh — ln[l + Fi('u)] 

4 J -OO ^ 



(50) 



(51) 



We have seen that the different states of the model are characterised by the zeroes of the 
Y-system and the cTa signs. To write down TBA integral equations for a given state one needs 
to know the structure of zeroes of the Y-system elements and the aa signs for the required 
state. These data can be read off (at least for large enough /) from the infinite-volume solution 
of the Y-system. In the SG model there exists a non-linear integral equation (DdV equation) 
which is used to describe the Bethe Ansatz states and calculate their energy and momentum 
[T3| IT^ IT5] . This equation contains only one (complex) unknown function and can easily be 
solved for large / in leading order (with exponential precision). The further advantage of this 
equation is that only those zeroes occur in it which give contribution to the energy (the set 
ri). The Y-system elements can be expressed by the unknown function of the DdV equation 
and this way the infinite-volume solutions can be obtained from the leading order solution of 
the DdV equation ^S]- It turns out that there are no 'complex holes' (C = 0) and the / — > oo 
solutions of the Y-system can be written in the |Im^| < 1 strip with exponential precision as 



-icosh f-^ 



(52) 



a= 2,...,p-l Dp+i 
YaiO = VaiO, { a= 2,...,p (53) 

a= 2, ...,p-2 Ap_2 
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where ^ ^ 

A(e) = (-1)^ i f[ e^>^(«-'^-+') + f[ e-'^^^-''--''^ \ 5 G {0, 1}, (54) 

la=l a=l ) 

the function 772(0 is defined by 

\{i + t)X{i-t) = l + m{0 (55) 
and the functions r]k{C) satisfy the Y-system equations 

%(e + - = [1 + %+i(0] [1 + %-i(0] (56) 

for = 2, 3, . . . and this determines^ ^fc(0 fo^' k > 2. 

It is possible to find the solution of ()56p explicitly. We note first that there is a class of 
solutions depending on a parameter q and a function B{^). Using these input data we first 
define 

k-l 

to(0 = 0, tfc(0 = $^g^i?[e + *(^-l-2j)] k = l,2... (57) 
and then it is easy to show that 

solve It is also true that 

The actual solution entering (58) for H soliton/fermion states corresponds to the choice 

H 

g = (-l)^ 5(^) = JJsinh^(e-M. (60) 

In the -Dp+i case the / 00 Y-system can be closed by defining 

Y,{0 = YMO = <0 = ^ff^y (61) 

The value of the parameter S in ()54p depends on whether one is considering the Sine-Gordon 
model soliton states or the masssive Thirring model fermion states as follows. 

SG : 6 = H (mod 2), MT : 5 = 0. (62) 



^?7i(^) = by definition. 
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From this it follows that the SG and MT models are the same when H is even and different 
when H is odd and this difference disappears in the infinite- volume limit. In the following 
we summarize the infinite-volume solutions of the Y-system for the ground state and the first 
excited state. 



H = Q\ (ground state) 



In this case we have a -Dp+i system for p > 2. From and we get 

A(0 = 2(-l)^ (63) 
Only the choice 5 = is physical. Further 

i?(0 = l, h{0 = k, Vk{0 = k''-l, KiO=p-l. (64) 



H = 1 ( one-particle state) 
Here, according to (fH^ we have two choices: 

6 = 1 for SG, (5 = for MT (65) 

and we have an A2p_i system for p >2 even and an v4p_2 system for p > 3 odd. The position of 
the hole is /ii = for the lowest lying excited state. This belongs to the the SG case and here 
we restrict our attention to this case for simplicity. From and (jSH) we find in this case 

A(0 = - {e'^(«+') + e-'''^^-'^ } (66) 

and from 



TT 



i?(0=sinh-(e) (67) 



we have 



cos"''' 



^ sinh ^ k odd 



^^(0 = <! lii ; (68) 

From this one can see that Va(0 has no zeroes for a odd and has a double zero at ^ = for 
even a values. All Ya{C,) are negative, except 1^(0 iii the A2P-1 case. With these analitical 
properties we have the following TEA equation for the first excited state. 

YaiO = i-lY^'"^ e-^-'™^^f« [r(0]'^^"'^" exp {/3,(0} a = l,...,p. (69) 

(a = 1, ... ,p — 2 for Ap_2-) ^a(0 /?a(0 even functions and Q;a(0 odd. It follows 
that the quantization conditions are automatically satisfied. Finally from (jKn|l and (jKT|) 

we get 

P« = 0, E« = M-^ c/Mcosh^ ln[l + ri(M)]. (70) 
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5 Liischer's formula 



In this section we make the assumption that there exist TBA equations describing the finite 
volume dependence of the 1-particle states in the SS- and sausage-models and that these are 
similar to the corresponding ground state TBA equations. For large volume the solution of the 
TBA equations can be obtained from Liischer's asymptotic formula, which is our starting point 
here. 

If a stable particle in a quantum field theory is enclosed in a box, its mass changes from its 
infinite-volume value due to the finite-size dependence of its self-energy. Liischer P derived a 
formula which describes the leading large- volume corrections to the mass of the lightest particle 
of the model in terms of the scattering amplitudes of the theory when periodic boundary 
conditions are imposed. This formula exists in all dimensions although it is really useful in two 
dimensional integrable models where the scattering amplitudes are exactly known. We discuss 
this formula in the simpliest two dimensional case, although it exists in higher dimensions 
too. If there is only one mass scale in the theory and there are no bound states, the leading 
large- volume correction of the mass scale is 

oo 

m{L) - M = -— / de cosh^ e-*'^™^'^'^a(^), (71) 
27r J 

— oo 

where 

^a{9) = E [-1 + ^atid + = + ^'^^^ + 'i^' ^^^^ 

b 

where n is the number of particles in the theory, m(L) is the mass gap in the theory enclosed 
in a box of size L with periodic boundary conditions, M is the infinite-volume mass and 

Qa{e)=Y.SM. (73) 

b 

When a quantum field theory is at finite temperature the virial coefficients of the pressure in 
the low temperature regime can be expressed by the scattering data alone [23]. The leading 
low temperature expression of the pressure is of the form 

oo 

p{T) = —nM I d9 cosh 9 6-^'°'^'^. (74) 
2tt J 

—oo 

Using the "modular transformation" j2] the pressure can be related to the ground state energy 
of the model in a box of size L with periodic boundary conditions: 

oo 

= -Lp (^i^ = j de cosh^ e-MLcoshe^ (75) 
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Using the results of (f73j) and (f7T|) the ground state energy {H = 0) and the first excited state 
energy [H = 1) of the model can be written in leading order for large L as 



oo 



E(^) = HM - Y / ducosh^u e-^^^^^"'^^!^ y^{u) H E {0,1} , (76) 

— oo 

where for the ground state energy 

H = 0, yi{u)=n, (77) 

and for the one-particle state energy 

H=l, y^(^u) = q,(^^iu + t)y (78) 

Let us assume that the TBA equations of our integrable model can be encoded in a Dynkin- 
diagram with one massive node and that the TBA equations of the first excited state of the 
model also exist. Then the energy of the ground state {H = 0) and the first excited state 
{H = 1) can be expressed in terms of the Y-system element associated to the massive node: 

oo 

E^"^ = HM - ^ j ducosh^u\n[l + Yi{u)] He {0,1}, (79) 

— oo 

where we assumed that the massive node of the Dynkin-diagram is indexed by one. Comparing 
(|7n|) with (f7n|) we get for Yi{u) in leading order 

Yi (u) = e-^ ^'^^ ^" Z/i (u) l = ML. (80) 

Applying this formula to the ground state {H = 0) of the SG-, SS- and sausage-models one 
gets the same result for Yi{u) in leading order as that coming from the large I solution of the 
corresponding TBA equations, (jHII), and respectively. Thus for the ground state the 
leading order coefficient of the "massive" Y-system element is determined by the leading order 
virial coefficient. 

In the previous section we have seen that (almost) the same Y-system describes all excited 
states of the SG model. (Small modifications occur at the end of the diagram in the odd charge 
sector of the model.) Calculating (f75jl using the scattering data (j21ISI) we get for the SG model 

q{e) = A{e) + B{e) = -e<'^') - e^>^(2*-v), (81) 

and from this we get 

yi(^) = - {e^xiiM) ^ g-ix(?-0} . (82) 
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This is the same as (jUBj) . which was obtained from the Bethe Ansatz solution of the modeL 
This agreement makes us confident that the method works if the assumption of the existence 
of the first excited state TBA equations is true. 

In the following we will assume that the TBA equations of the first excited state of the SS- 
and sausage-models also exist and they are described by almost the same Y-system as for the 
ground state. (The position of the massive node is the same and small modifications can occur, 
like in the SG model, at the end of the diagram.) 

Making this assumption the leading order expression of Yi{u) can be read off (|78|) using the 
scattering data of section 2. We get for the leading order expression for Yi{u) in the SS-model 

yiio = -m (83) 

where 

= gix«+i) _e-ix(«-i)^ (84) 



5^^) = _e-ix«-i). (85) 

In the sausage-model one gets different results for the charged and for the neutral particle 
states. The leading order result for the charged particles is 

where 

BiO = sinh + ^) ■ si^h ^(e - (87) 
and for the neutral particle states we have 

yoiO = ^ _ {35(0 - sinvrA ■ sin27rA} ^ y^iO- (88) 

Starting from these large / expressions of the "massive" Y-system elements, the / —>■ oo solution 
of the full Y-system can be obtained using the Y-system equations pOjl recursively and the 
structure of the zeroes of the Y-system elements can be determined. 



6 The / oo solution of the one particle Y-systems 

In this section we will calculate the / oo solution of the Y-systems for the first excited states 
in the SS- and sausage- mo dels. 
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The SS-model 

Starting from the massive node expression (jH!?|) and applying the Y-system equations (jHn|l 
recursively the infinite-volume solution of the Y-system can be obtained: 

n(0 = -e-'^°^'^«&(0-KO, (89) 

Yk{0=Vk{0 k = 2,...,p, (90) 

Yk{0 = Vk{0 k = 2,...,p. (91) 

Here rjk{^) is of the form with q = I and 

i?(0 = sinh^, (92) 
2p 



nk 



sm ^ 

sin ^ 2p 

and 772(0 is determined by 



4(0 = —^ sinh^ (93) 



K^ + OK^-0 = -[1 + ^2(0]- (94) 

Using the reflection symmetry 

tp+niO = tp-niO, Vp+niO = Vp-n{0 (95) 

the diagram is closed by the relation 

VpiC + ^)Vpi^-^) = [l + Vp-iiO?- (96) 
The same expressions ()92ll9(i|l hold for the fjk{^) quantities with the replacements 

This Y-system is encoded by the diagram of Figure 5. Modifications with respect to the 
ground state problem occur at both ends of the diagram, similarly to the SG case. 
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The sausage-model 

Starting from the massive node expression (jHH|) and applying the Y-system equations (jHn|l 
recursively the infinite-volume solution of the Y-system of the charged particle states can be 
obtained: 

yi(0 = -e-''=°^'^5«r^2(0, (9T) 
Yk{0=Vk{0, k = 2,...,N-l (98) 

Y,,iO=yN+iiO = <0, (99) 
where rik{^) and k,{^) can be written in the form of (jKHj l?)T|l with 

q=l r/2(0=y+(0, (100) 

sinAvr/c , , /. ^ , 

ikiO = 7r~- T -^^^ ~ o (101) 

2 sm ttA 2 

"^^^ " ~cosh2^(e + 0-cosh2^(e-«)' 
The corresponding TBA diagram is exactly the same as for the ground state (Figure 3). 

In the next two sections we will use these / oo solutions of the excited state Y-systems to 
determine the structure of zeroes of the Y-system elements, which is necessary for transforming 
the Y-systems into TBA integral equations. The / oo solutions will also be used as starting 
functions in the iterative numerical solution of the TBA integral equations. 

7 One-particle TBA equations of the SS-model 

In this section the one-particle TBA equations of the SS-model will be written down by using 
the analitical properties of the Y-sytem elements determined by the infinite-volume solution 
()89ll96p . From this we see that all Ya (^ ) functions have a double zero at ^ = and have no other 
zeroes and that all Ya{oo) > 0. From these properties using one can derive the following 
integral equation for the one-particle state. 

YaiO = e-'^-^°^^5? e^"^^^ « e {l,2,..,p, 2,..,p}. (103) 

The "quantization conditions" Fa(±?) = —1 are satisfied automatically because Ya{i) and /9a(0 
are even functions and aa{Cj are odd. 
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p p-1 p-2 2 1 2 p-2 p-l P 



Figure 5: Dynkin-diagram associated witli tlie SS-model Y-system (excited 
state) . 



12 3 4 

Figure 6: Dynkin-diagram associated with the 0(4) model Y-system. 



The 0(4) model is the p,p ^ oo Umit of the SS-modeL In this hmit the Y-system consists 
of infinitely many components and the TBA diagram becomes symmetric to the massive node: 

n(e) = n(0 A; = 2,3,... (104) 

The corresponding infinite TBA diagram is depicted in Figure 6, where the oriented double 
line at the beginning of the diagram means 

/i2 = 2, hi = 1. (105) 

The infinite- volume solution of the 0(4) model Y-system is of the form 

n(0 = ^^e^ A. = 2,3,... (107) 

8 One-particle TBA equations of the sausage-model 

In this section the one-particle TBA equations of the sausage-model will be written down by 
reading off the analitical properties of the Y-sytem elements of the infinite-volume solution 
()97II102|) . We will consider only the generic N > 5 case. From ()97H102j) the signs at infinity are 

rikioo)>0 k = 2,3,...,N -2, (108) 



rjN-i{oo) < 0, 



K,(oo) < 0. 



(109) 
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The zeroes of the infinite-volume Y-system are solutions of the equation 



V2 : 


: ±H3, 


m ■■ 


: 0,0; ±H,, 




: ±Hs^i; ±H. 




: ±if7v-2, 


K : 


: ii^AT-i. 



, , k sin2A7r 

coshA7ri7fc = — --— , 110 

/ Sm KATT 



H2 = 0;H2<H3< ... < Hn^i. (Ill) 
All HkS are real and the zeroes of the infinite- volume solutions of the Y-system are 

(112) 
(113) 

s = 4,..,iV-2 (114) 

(115) 
(116) 

Using the above analitical properties one can derive the following one-particle TBA equations 
for the sausage-model. 

Y2{0=r{i-H,)r{i + H,)e^-^^\ (117) 
n(0 = e-'^°^"^« 1^2(0, (118) 

n(0 = r(C - if,_i) r(e + i/.-i) r(e - iJ.+i) r{i + H,+^) eM^\ s = 3,..,N-2, (119) 

Yn-i{0 = - Hn-2) r(e + Hn^2) e^^-^^^\ (120) 

Yn{0 = yN+iiO = -r{^ - Hn-i) r(e + H^-i) e^-(«). (121) 

These equations must be supplemented by the quantization conditions (j44|) . The first two 
quantization conditions Yi{±i) = ¥2(^1) = —1 are satisfied automatically due to the fact that 
H2 = exactly. The rest is of the form 

^{H, - H,^,) + 7(/7, + Hs-i) + i{Hs - H,+i) + i{H, + Hs+i) - a,{Hs) = 27rM„ (122) 

s = 3,...,A^-2 

-i{HN^i - Hn-2) + i{Hn-i + Hn-2) - aN-i{HN-i) = 2tiMn-i, (123) 
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where 

Ms : half-integers s = 3, . . . , iV - 1 (124) 

and can be determined from the infinite- volume solutions ()97lll02j) . 

The 0(3) model is the N oo limit of the sausage-model, where the TBA diagram becomes 
infinite. In this case the infinite-volume solution becomes 

Fi(e)-e-'™^^i«^^, (125) 



n(o = ''i, , ^ = 2,3,... (126) 



( p - 1) ie - Dk-i) ie - Dk+ i) 
[e + {k- 1)2] ■ [e + ik + 1)2 

where 



3 

and the zeroes of the infinite-volume Y-system are of the form 



A;2 - 4 

Dk = (127) 



Hu = ^/dI k = 2,3,... (128) 



9 Numerical results 

Our aim in this paper is to calculate the finite volume mass gap for the 0(3) and 0(4) nonlinear 
cr-models. Since these are the N oo limit of the sausage-model and the p,p —>■ oo limit of 
the SS-model respectively, one has to solve the TBA integral equations introduced in sections 
3, 7 and 8. 

Concretely, for the ground state energy {H = 0) we have to solve of section 3. For 
the first excited states {H = 1) in the SS-model we have to solve p()3|) of section 7. Finally 
the first excited state problem in the sausage-model requires the solution of the TBA integral 
equations pi7lll2l|l together with the quantization conditions p22lll24|) . After having solved 
the integral equations (|79|l can be used to calculate the energies and E(i). 

For finite (sausage-model) and finite p,p (SS-model) the TBA problem can easily be 
solved numerically by iteration. For H = and also the H = 1 case of the SS-model this is 
completely straightforward. As usual, the / oo solution can be used as starting point for 
the iteration and the procedure converges rather rapidly. The excited state problem for the 
sausage-model is more involved^ since here one step in the iteration includes the calculation 
of the integrals occuring in ()117H12l|) together with the calculation of the roots of ()122H124|) . 
•^This is similar to the H = 2 problem in the SG case. 
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A 
1 v 




u 


6 




U{6) model H = {) 


140 


140 


0.1 


10 /lO 




C^Jj model n = l 


oO 


oO 


0.05 


iU /iU 


1 n— 7 /I n— 8 
10 /lO 


U[4:) model H = [) 


1 /I n 
i4U 


1 /I n 
i4U 


O.i 


1 n— 6 / 1 n— 7 
lU /lU 




0(4) model H = 1 


100 


60 


0.1 


10-8 





Table 1: Parameters for numerical calculation of the ground state and first excited state energies in the 
0(3) and 0(4) nonlinear (T-niodels. Where two values are given for e and ei the first one refers to the range 
0.001 < I < 0.1 and the second one to the range I > 0.1. 

Again, the starting point of the iteration procedure is given by the / — » oo solution, both for 
the Y-system functions and the position of the zeroes Hg. Before the iteration is started, the 
half-integers in ()124j) have to be calculated using the large volume solution. We found that in 
all cases all half-integers are equal to one half. 

To calculate the a-model limit one has to take large N (or p, p) values and extrapolate. We 
adopted a slightly different approach: we studied a cutoff oo system. The cutoff system for the 
first excited states (with cutoff u) is obtained by considering the N = oo {p = p = oo) limit of 
the TBA problem but "freezing" Ya{C,) for a > u at the large volume limit solution, which is 
given by ()126|) and ()107|) for the 0(3) and 0(4) model respectively. The cutoff infinite system 
for the ground state problem is defined analogously. In this way we found faster convergence, 
probably because for the cutoff infinite system the starting point of the iteration (also for a < u) 
is the (T-model limit and thus closer to the final solution. 

In Table 1 we summarized the values of the parameters we used for numerical determination 
of the ground state and first excited state energies for the 0(3) and 0(4) nonlinear a- models. 
Here A is the rapidity cutoff, 6 is the length of the intervals used in Simpson's formula, e is 
the relative precision of the numerical results for the energies and finally ei is the accuracy of 
the solution of the quantization conditions p22j) . To the required precision (e) there is no need 
to take better A, z/ or 5 values^ than those given in Table 1. With these parameter values the 
required precision is achieved after a few thousand iterations, which are typically completed in 
less than an hour on a PC. The only exception is the excited state problem for the 0(3) model, 
which requires a few hours of CPU time on a PC. Our numerical results are summarized in 
Tables 2 and 3. All energies are given in units of the infinite volume mass M. 
''At least for the range I > 0.001 we consider here. 
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I 






0.001 


-913.954(1) 


-406.23(1) 


0.003 


-299.5016(3) 


-111.903(3) 


0.01 


-87.6357(1) 


-23.643(1) 


0.03 


-28.27949(3) 


-3.8325(3) 


0.1 


-8.006985(1) 


0.77718(1) 


0.3 


-2.3890980(3) 


1.235363(3) 


1.0 


-0.4862496(1) 


1.084208(1) 



Table 2: Numerical results for ground state and first excited state energies in the 0(3) nonlinear cr-model. 



I 






0.001 


-1343.408(1) 


-901.28159(1) 


0.003 


-438.1506(3) 


-272.740308(3) 


0.01 


-127.2263(1) 


-69.838028(1) 


0.03 


-40.60919(3) 


-18.2468766(3) 


0.1 


-11.273364(1) 


-3.0041089(1) 


1/8 


-8.8346989(8) 


-1.91603183(8) 


1/4 


-4.0487295(4) 


-0.00446984(4) 


1/2 


-1.7404694(2) 


0.71072801(2) 


1 


-0.6437746(1) 


0.93839706(1) 


2 


-0.16202897(1) 


0.99233406(1) 


4 


-0.015(32571(1) 


0.99965327(1) 



Table 3: Numerical results for ground state and first excited state energies in the 0(4) nonlinear cr-model. 
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10 Monte Carlo and perturbative results 



We are now in a position to be able to compare our numerical results with those of Monte 
Carlo simulations and perturbative calculations. Checking our results using asymptotically 
free perturbation theory (PT) in the small volume {I — > 0) limit is especially important since 
our construction is based on Liischer's formula and the large volume (/ — > oo) solution. 
Perturbative calculations for the finite volume mass gap of the 0{n) nonlinear cr-model 

z{l) = LM{L) = I [e^^^ (l) - (/)] (129) 

are available up to three loop order [2E]- The results are best presented in the form of the 
asymptotic expansion 

7r(l + A) f A M3 1 , , 

where A = l/(n — 2) and the (inverse) running coupling x is the solution of the equation 

x-Alna; = ln(— ^ ) . (131) 



The perturbative lambda parameter used here is the Finite Volume lambda, which is best suited 
to the problem j2E] and is related to the conventional by 



AFv = ^AMs = ^(8j r(l + A)M. (132) 

Here 7 is Euler's constant and in the second equality we have used the exact value of the 
M/Ay^ ratio, which is available for this family of models fI7\ . 
The coefficient of the three-loop contribution is [2E] 

«3 = f + fA + |A'H-(|i-i)A3, (133) 

where 

Xi = -1.2020569, X3 = 23.6, 

(134) 

X2 = -3.63, X4 = -5.2123414. 

In terms of the physical volume / we have to solve 

X — Alnx = In — ui, (135) 

where 

=7-ln47r + Aln(^0 +lnr(l + A). (136) 
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Figure 7: Finite volume mass gap of the 0(3) model. Comparison of numerical solution of the TBA integral 
equations (dots) to three- loop perturbation theory (solid line). The two- loop perturbative curve (dashed line) 
is also shown. 

The comparison of our numerical results to the three-loop perturbative predictions is shown 
in Figures [3 and |H1 It is reassuring to see that our results agree very well with asymptoically 
free PT in the small volume regime. 

In Table 4 we collected all available data on the finite volume mass gap for the 0(4) model. 
In addition to our numerical results the three- loop perturbative results (for small volumes), 
results of MC simulations j2H] and (for large volumes) the values corresponding to Liischer's 
formula are given. 

Some MC results are available also for the 0(3) model. We intend to discuss how they 
compare to the results of the TBA calculations presented in this paper in a future publication. 
Here we only mention that we have calculated the value of the "step scaling function" [3] 
(j(2,uo) at the "canonical" point Uq = 1.0595. We found 

a(2, Mo) = 1-261208(1), (137) 

which seems to support the results of the form-independent fit in P^ . 
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t 


J. UjV 


PT 


A/TP 


1 nO/^rlOT* 


U.UUl 


U.44212o 


A /I /I OA A 'V 

0.442097 






(J.UUo 


A A ACOO 1 


A /I A/^ 1 'vr' 

0.49oi7o 






A A1 


A C 'TO O O O 

U. 573883 


0. 573700 






A AO 


U.d7U8d9 


O.o70o0o 






A 1 

U.l 


A OOi^ AO !X 

(J.8zd9z55 


A OO/^ 1 OA 

0.82oi30 




A A TOO 

U.4733 


1 /o 

1/8 


U. 8648334 


A O O O O 1 

0.8o3821 


A OGO /0\ 

U. 863(2) 


A CT O/^ C 

0.5265 


1/4 


1 A1 1 A/^ /I A 

1.U11Ud49 


1 A A O n O 

1. 00868 


1 A1 1 /0\ 

l.Ull(2j 


A '"TOGO 

0.7368 


1/2 


1 OO C C AO '7 

I.zz55y87 


1 O 1 O/^ '7 


1.228(2j 


1 AO AC 

1.0395 


1 


"1 r o o "1 "1 




i.bo 1( 1 j 


1. 19 li 




z.otJo / 2b i 




2.3(Jy(i(Jj 


2.2909 


1 


1.0611160 




1.132(10) 


1.0(307 



Table 4: Results for the finite volume mass gap z{l) in the 0(4) nonlinear tr-model. 

11 Summary and conclusion 

In this paper we proposed TBA integral equations for the 1-particle states in the sausage- and 
SS-modcls and for their (T-model limits, for the 0(3) and 0(4) nonlinear cr-models respectively. 
The excited state TBA systems are based on analogy with the corresponding problem in the 
Sine-Gordon model and the solution in the large volume limit, which can be obtained from 
Liischer's asymptotic formula. Combining the 1-particle TBA systems with those corresponding 
to the ground state we can calculate the exact value of the mass gap numerically for the sausage- 
and SS-models, and, by extrapolation, also for their cr-model limits, which correspond to infinite 
TBA systems. We have proposed a somewhat different, more efficient method to treat these 
infinite TBA systems: instead of taking larger and larger sausage-model (or SS-model) TBA 
systems and extrapolating, we consider directly the infinite system, with a cutoff that removes 
the remote TBA nodes. Considering the cutoff infinite system instead of the original problem 
leads to faster convergence of the iteration procedure and produces numerically precise results 
for the cr-model mass gap already at moderate values of the cutoff parameter. 

Having computed the mass gap for the 0(3) and 0(4) models we can compare the results 
to those of lattice Monte Carlo simulations and perturbation theory. We have observed perfect 
agreement taking account of all available data. Since the cr-models are asymptotically free, the 
perturbative results are reliable at small volumes. On the other hand our TBA systems are 
based on Liischer's large volume asymptotic formula and hence the very good agreement of 
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Figure 8: Finite volume mass gap of the 0(4) model. Comparison of numerical solution of the TBA integral 
equations (dots) to three- loop perturbation theory (solid line). The two- loop perturbative result is too close to 
the three-loop result and is not shown here. 

the perturbative results with our numbers for small volumes indicate that all our assumptions 
leading to the excited state TBA systems are vahd. The lattice Monte Carlo data are for the 
intermediate volume range and the results agree with our mass gap values within the Monte 
Carlo errors. Thus we have numerically checked the correctness of the proposed TBA equations 
for a large volume range between I ~ 1 and I ~ 10~^. 
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